A limiting strategy for the back and forth error compensation and correction method for solving advection equations

We further study the properties of the back and forth error compensation and correction (BFECC) method for advection equations such as those related to the level set method and for solving Hamilton-Jacobi equations on unstructured meshes. In particular, we develop a new limiting strategy which requires another backward advection in time so that overshoots/undershoots on the new time level get exposed when they are transformed back to compare with the solution on the old time level. This new technique is very simple to implement even for unstructured meshes and is able to eliminate artifacts induced by jump discontinuities in derivatives of the solution as well as by jump discontinuities in the solution itself (even if the solution has large gradients in the vicinities of a jump). Typically, a formal second order method method for solving a time dependent Hamilton-Jacobi equation requires quadratic interpolation in space. A BFECC method on the other hand only requires linear interpolation in each step, thus is local and easy to implement even for unstructured meshes.